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Abstract 

Given a set of unit-disks in the plane with union area A, what fraction of A can be covered by 
selecting a pairwise disjoint subset of the disks? Rado conjectured 1/4 and proved 1/4.41. Motivated by 
the problem of channel-assignment for wireless access points, in which use of 3 channels is a standard 
practice, we consider a variant where the selected subset of disks must be 3-colourable with disks of the 
same colour pairwise-disjoint. For this variant of the problem, we conjecture that it is always possible to 
cover at least 1/1.41 of the union area and prove 1/2.09. We also provide an O(n^) algorithm to select 
a subset achieving a 1/2.77 bound. 

1 Introduction 

Richard Rado studied the following problem: What is the largest ci such that, given any arrangement of 
unit-disks D in the plane, we can always select a pairwise disjoint subset of disks that cover at least a 
fraction ci of the area of the union of D? Clearly ci < |, corresponding to the case shown in Figure [T] (left) 
where a large number of unit-disks share a very small intersection — the common intersection prevents us 
from selecting more than a single disk, while the union area of all disks approaches 4tt. Rado proved a lower 
bound of ci ^ ~ [H] and conjectured the lower bound ci ^ j. 





Figure 1: Given a set of disks arranged in a circle with a very small mutual intersection, the most area we 
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can cover is 4^ = i with one colour (left) and approximately with 3-colours (right). 



In this paper, motivated by the problem of channel-assignment for wireless access points, in which use 
of 3 channels is a standard practice, we focus on a variant of the above problem in which the disjointness 
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constraint on the selected subset of disks is relaxed slightly. Given an arrangement of unit-disks D in the 
plane, we want to find the largest C3 such that we can always select and 3-colour a subset of the disks C 
such that their union area covers at least a fraction C3 of the union area of under the constraint that 
same-coloured disks must be pairwise disjoint. We will prove that, for any given arrangement of unit-disks, 
C3 ^ 2l)9' This result is stated formally in Theorem [T] 

Theorem 1. Let D he a collection of unit-disks in the plane with union area A. For C a 3-coloured subset 
of D with same- coloured disks pairwise disjoint, let Ac denote C's union area. We can always select a C 
such that ^ ^ 2^ '^'^'^ thus C3 ^ . 

Note that C3 < : corresponding to selecting three minimally pairwise intersecting disks in an arrange- 
ment as shown in Figure [1] (right). We conjecture that this bound is achievable for any arrangement of 
unit-disks. 

Before presenting our proof of Theorem [1] we will review Rado's proof and prove the weaker result that 
C3 ^ Yff 1 formally stated in Theorem O We include this proof as a stepping stone to Theorem [1] and 
because it leads to an 0{n'^) algorithm for the disk selection problem. 

Theorem 2. Let D be a collection of unit-disks in the plane with union area A. For C a 3-coloured subset 
of D with same-coloured disks pairwise disjoint, let Ac denote C's union area. We can always select a C 
such that ^ ^ • 

The rest of this paper is organized as follows. In Section [2] we present our motivation for exploring this 
problem and discuss some related problems. In Section [3] we review Rado's proof and prove Theorem [21 In 
Section m we prove Theorem [1] In Section [5] we present an O(n^) time algorithm to select a subset which 
achieves our 1/2.77 bound from Theorem [5] and prove that the lattice-positioning step in our algorithm is 
3SUM-hard. Finally, Section [6] discusses bounds for fc-colours when fc 7^ 3. 



2 Motivation and Related Work 

In this section we describe a wireless network deployment problem that motivates our 3-coloured disk packing 
problem, and discuss other disk packing results related to ours. 

Wi-Fi (IEEE 802.11) wireless networks are becoming a ubiquitous feature in modern businesses, univer- 
sities, parks, etc. In a typical Wi-Fi deployment scenario, a set of candidate locations are determined for 
wireless access points ( APs) . A subset of the candidate locations must be chosen along with a channel assign- 
ment for each installed AP in order to maximize the area covered by the wireless network while minimizing 
interference. Interference occurs when two APs using the same channel are within range of one another, 
preventing users in range of both APs from communicating with either. For a more in depth discussion 
of wireless network deployment, see the introduction to [3]. We make the simplifying assumption that the 
coverage area of each AP is a disk, though in practice coverage areas may be more irregular. We also assume 
that every AP uses the same power level — i.e. that all the disks have the same radius. This is not entirely 
unjustified, as the problem of power control for Wi-Fi network devices is a complex research area in its own 
right pj. Finally, the set of potential locations for APs is limited to those that are physically possible and 
aesthetically acceptable. We model this as a finite set of disk center points. Under these conditions, the 
wireless network deployment problem becomes a fc-coloured disk packing problem, where k is the number of 
allowable channels. To justify k = 3, we will discuss the constraints on channel selection in Wi-Fi networks. 

The IEEE 802.11 standard establishes a number of requirements on the radio frequency characteristics of 
802.11 hardware. The 2.4-Ghz band used by 802.11 devices is broken into 11 channels for the North American 
domain and 13 channels for the European domain. Unfortunately, the number of effective channels is much 
lower due to interference between channels with center frequencies close to one another. A study by wireless 
hardware maker Cisco Systems recommends that wireless network deployments only use three channels (1, 
6, and 11 for the North American domain), finding that even a four-channel scheme can cause unacceptable 
degradation of service in systems with a high volume of users [10] . 
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A more sophisticated formalization of the deployment problem allows disks assigned to the same channel 
to overlap but only counts the area where there is no interference — i.e. the area of the set of points covered 
by only one disk on some channel. In terms of colouring, the problem is to colour a subset of the given disks 
to maximize the area of £ | for some colour, point p is in exactly one disk of that colour}. We call this 
the 1-covered area. Asano et al. [5] proved that it is always possible to achieve approximately 1-covered 
area using only one colour. This model has also been considered with respect to two other optimization 
problems. For the problem of maximizing the 1-covered area using one colour, previous work has focused 
on approximation algorithms (though no proof yet exists, it is suspected that this problem is NP-hard). In 
[5], Asano et al. present a 5.83-approximation algorithm with polynomial runtime. In [9], Chen et al. show 
that the problem admits a polynomial time approximation scheme when the radius of the largest disk over 
the radius of the smallest disk is a constant. 

Another well-explored optimization problem is conflict-free colouring — here the goal is to minimize the 
number of colours needed to 1-cover the whole area, i.e. the union of the given disks. Even et al. \lT\ prove 
that 0(log n) colours are always sufficient and sometimes necessary for any given disks of general radii. Alon 
et al. [1] have shown that, if each disk intersects at most k others, then 0{log'^ k) colours are sufficient 
for a conflict-free colouring, improving the bound from |11) when k is much smaller than n. There is also 
work on online algorithms for conflict-free colouring |12j , and on conflict-free colouring of regions other than 
disks tl4j. 

3 The Basic Bound 
3.1 Rado's Proof 

The idea behind our proof of the basic bound is similar to that used by Rado in [16]. In Rado's proof, a 
regular triangular lattice of side length 4 is positioned over the set of disks D and for each lattice point that 
falls in UD one disk containing that lattice point is selected. The side length of the lattice guarantees that 
disks selected in this manner will be pairwise disjoint (see Figure [21). Thus, supposing we can prove a lower 
bound of k on the number of lattice points intersecting any given set of disks, we immediately obtain a lower 
bound of ci ^ ^ . We get such a lower bound by applying Lemma [3] which states that we can position the 
lattice to contain at least ^ points in UD, where a is the area of a triangle in the lattice. Since a = 4-\/3 
for our lattice, we can position the lattice to contain at least k — J^p- points in UZ) and therefore ci ^ ^ . 




Figure 2: A triangular lattice ensures that selected disks are pairwise disjoint. 



Lemma 3 ^16j . Given a region of the plane G with area A, and a triangular lattice in which each triangle 
has area a, the lattice can always he positioned such that it contains ^ points in G. 

Proof. Rado's proof [16] uses the concept of the "fundamental cell" of a lattice — for a regular triangular 
lattice the fundamental cell F consists of a pair of adjacent triangles (see Figure [3]). Given an arbitrary 
placement of the lattice, each triangle of the lattice can be translated to F along with whatever part of 
G is contained in the triangle. The translated parts of G may "overlap" in F — there may be points on 
the fundamental cell intersecting multiple translated portions of G. Supposing a point p in F intersects k 
translated portions of G, then repositioning the lattice such that p is a lattice point ensures that k lattice 
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points intersect G. In this case, we refer to k as the depth of point p in the fundamental cell. Since the area 
of G is A, the total area of all portions of G translated to F is clearly also A. Therefore, we have portions 
of G with total area A translated to a region of area 2a and a point of depth at least ^ must exist. □ 

3.2 A Bound of cg > 1/2.77 

Theorem 2. Let D be a collection of unit-disks in the plane with union area A. For C a 3-coloured subset 
of D with same- coloured disks pairwise disjoint, let Ac denote C's union area. We can always select a C 
such that ^ ^ 277- 

Proof. To solve our variation of the problem we use a finer triangular lattice with side length ^-j^- The 
points of the lattice are 3-coloured such that no two lattice points of the same colour are adjacent. For any 
placement of the lattice, select a subset C of Z) as follows: for each lattice point p in the union of D, select 
a disk containing p and assign the disk the colour of p. The side length of the lattice ensures that no disk 
contains two lattice points so the selection and colouring are well-defined. It also ensures that disks assigned 
the same colour are pairwise disjoint (see Figure [3]). 




Figure 3: A finer 3-coloured lattice ensures only that same-coloured selected disks are pairwise disjoint. 



Also observe that, by Lemma [31 we can position the lattice to intersect the union area of D in at least 
points, so |C| ^ ^^8^- While same-coloured disks in C are pairwise disjoint, differently coloured disks 
may not be, so |C|7r is only an upper bound on Ac- 

To derive a lower bound we will partition the union of C using the triangular lattice's Voronoi tessellation 
which has regular hexagonal cells of side length | and vertices at the barycenters of the triangular lattice 
(see Figure E]). Suppose disk d G C contains lattice point p which lies in hexagonal cell h. If we count only 
the area of d n ft., and sum over all d, this gives a lower bound on Ac- Thus if we establish a lower bound 

A on the minimum possible area ofdnh then Ac > |C|A ^ ^A. In LemmaH which we will prove in 
Section 13.31 we show that A « 1.6645. From the lower bound on Ac we reach our desired lower bound on 
C3 of 

^Ac Vs. 1 

c-A > — — = A ^ 

A 8 2.77 

□ 



3.3 Minimum Disk-Hexagon Intersection 

Lemma 4. Given a regular hexagon h with center point X and side length |, and any unit-disk d containing 
point X , the minimum area of intersection A between h and d is approximately 1.6645, or more precisely 

. ^ , yiT , TT 1 ^ fbVs-Vn 

A = H H arctan 



36 12 2 2 \5 + VnV3 
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Figure 4: The Voronoi tessellation of the triangular lattice points forms a grid of regular hexagonal cells. 



Proof. We use elementary geometry to argue that the minimum area of intersection is achieved by a disk d 
with X on its boundary. Then, parameterizing by the angle 6 between the horizontal axis and the ray from 
X to the center of d, we use the symbolic geometry package Geometry Expressions to give a formula for the 
area of intersection and use Maple to compute O's of the first derivative, finding that the minimum is as stated 
above, and occurs in the configuration shown in Figure [S] Further details are included in Appendix \X[ □ 

Our proof of Lemma|4]also shows that the lower bound Ac > \C\A is tight, as shown by the example in 
Figure [5] where the union of C is exactly partitioned by the hexagons and each disk intersects its hexagon in 
the minimum area A. However, note that this does not mean that our bound on C3 is tight. In particular, 
for the example shown in Figure [5] we can capture the whole area by 3-colouring the disks. 




Figure 5: In this arrangement of selected disks, the lower bound on the contribution of each disk is realized. 



4 Deriving a Better Bound 

In our proof of Theorem [5] we only counted the minimum intersection of a disk with its selecting Voronoi cell 
(the Voronoi cell of the lattice point that selects the disk). Suppose we commit to this accounting scheme. 
We can clearly maximize the intersection of the disk with its Voronoi cell by centering the selecting lattice 
point in the disk. This suggests that we can use a weighting function that prefers placing a lattice point in 
the center of a disk, and do a more intricate analysis of the contribution of each disk. 

Looking at this from another direction, in the proof of Theorem [5] we optimized the number of disks 
selected, rather than the area of the intersection between selected disks and their selecting Voronoi cells. 
This approach can be improved because, among all subsets of disks that can be selected using a lattice, the 
largest subset of disks does not necessarily cover the largest area (e.g. see Figure [6] where selecting the three 
intersecting disks using the lattice shown gives a subset with union area 3A « 4.99, while an alternate lattice 
positioning that selects two disjoint disks gives a subset with union area 2it « 6.28). 

The above points suggest that we can improve our bound on C3 by using a more sophisticated criterion for 
lattice positioning based on the area contributed by selected disks rather than the number of disks selected. 
We will use this approach to prove Theorem [TJ 



5 



Figure 6: The lattice positioning shown selects the maximum number of disks, however an alternate lattice 
positioning can select fewer disks and cover more area. 



4.1 A Bound of C3 > 1/2.09 

Theorem 1. Let D be a collection of unit- disks in the plane with union area A. For C a 3- coloured subset 
of D with same- coloured disks pairwise disjoint, let Ac denote C's union area. We can always select a C 
such that ^ ^ 21)9 '^^'^ thus C3 ^ 2ijg ■ 

Proof. To prove this better bound, we use a weighted version of Lemma [31 Whereas Lemma [3] proved that 
a lattice L can be positioned to intersect \JD in at least a certain number of points, we want to prove that 
L can be positioned to intersect UD with at least a certain weight. 

We now define the weight function. For point p, let H{j)) be the regular hexagon of side length | 
centered at p. Thus H{p) is the Voronoi cell of p if our triangular lattice is translated to include point 
p. For p g U-D, let d{p) be the disk containing p whose intersection with H{p) has maximum area. Now 
let w{p) — area{H (p) D d{p)) for p € UD, w{p) = otherwise. Then w{p) measures the area contributed 
by including p in the lattice (given our method of choosing disks based on the lattice, and our accounting 
scheme of counting only the area of the disk in the Voronoi cell). 

We want to choose a lattice L to maximize W{L) = X^pGL ""^(p)- Proving a lower bound on the maximum 
value of W{L) for lattices of the type used in our proof of Theorem [2] will give us a lower bound on C3. 

As in the proof of Lemma [3] we consider translating portions of UD onto a fundamental cell. After 
translation, each point p on the fundamental cell is assigned the sum of the weights of p's intersection 
with each translated portion of UD. Supposing that a point p on the fundamental cell has weight to, then 
positioning L such that p is a lattice point ensures that W{L) = m. Thus, if we can prove that the total 
weight of the portions of UD translated to the fundamental cell is B, then a point in the fundamental cell 
(and therefore a lattice positioning) with weight at least must exist. 

We cannot compute exactly the value of B, which is the integral of w{p) over p G UD, but we can obtain 
a lower bound. Our first step is to replace the true weight function w{p) by a lower bound, wi{p), in which 
we replace the hexagon by its inscribed circle (see Figure [7]). Specifically, wi{p) is the area of the intersection 
of two discs: the disk of radius 2/V3 centered at p, and the unit-disk in D whose center is closest to p. Note 
that wi{p) depends only on the distance, r, from p to the nearest disk center in D. We will overload the 
notation and define wi{r), for r g [0, 1], to be the area of the intersection of a unit-disk and a disk of radius 
2/^/3 whose centers are distance r apart. We note that wi{r) is a continuous, decreasing function, and that 
we can write an explicit formula for it: 
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Figure 7: Considering the intersection between a selected disk and the largest disk that can be contained 
in the disk's Voronoi cell provides a simple lower bound for 'w{p) in terms of distance from p to the nearest 
disk center. 



ifO<r<^ 



Mr) = ' arccos (i ^) + | arccos U i:!±iM' 



^V("" + ' + V5)('- + l-75)('--l + 7l) (^ + l + 7l) ifTS-K'^^l 

We will express B as an integral in terms of variable r. Note that the points on the boundary of UD are 
precisely the points at distance 1 from the closest disk center. More generally, we can capture the points 
that are distance r from the closest disk center as follows. For unit-disk d G D, let dr be a disk of radius 
r at the same center. Let Dr be {dr : d G D}. Then the points that are distance r from the closest disk 
center are precisely the points on the boundary of UDr. Let p{r) be the length of the boundary of UDr- As 
w{r) > wilr) > 0, we obtain: 

1 1 

B = J p{r)w{r) dr > J plr)wi{r) dr 



We prove a lower bound on the latter integral that eliminates p{r) . 
Lemma 5. _B = p{r)w{r) dr > p{r)wi{r) dr > 2 A rwi{r) dr. 

After this we are done: we plug in the expression for wilr) and evaluate the integral using Maple to 
obtain a lower bound on B of 2.2Q1A. Therefore a lattice positioning with weight at least ^B ^ must 
exist, and we get our desired lower bound of C3 ^ 2l)9- 

It remains to prove the Lemma. 

Proof of Lemma\^ We want to prove Jq p{r)wi{r) dr > 2 A rwi{r) dr, equivalently: 



[p{r) - 2rA]wi{r) dr > 

Let A{r) be the area of UDr, which is related to the perimeter by the fact that p{r) — A'{r). Define 
/(r) = A{r) - r'^A. Recall that A = A{1) is the area of UD. Note that /(O) = /(I) = 0. Now /'(r) = 
A'{r) — 2rA = p{r) — 2rA, and the inequality we want to prove is transformed into f'{'r)wi{r) dr > 0. 

We apply integration by parts, noting that /' is continuous, and that w'l is continuous except at — 1. 
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f'{r)wi{r)dr = f{r)wi{r) 



1 1 
f{r)w'i{r)dr = — f {r)w[{r) dr 



We will prove below that A{r) > r^A. Thus /(r) > for all r £ [0, 1]. Now w;(r) is a decreasing function, 
so w'i{r) < for all r € [0, 1]. The integral of a negative function is negative, and this completes the proof. 

□ 

□ 





Figure 8: The union of a set of disks before (left) and after a radius scaling has been applied (right). 



Claim 1. Given a collection of unit-disks with union area A, if we scale the radius of each disk by r [0, 1] 
then the union area of the scaled disks will be at least r^A. 

Proof. We want to scale the radius of each disk by r. We accomplish this in two steps. First we scale the 
whole plane by a factor of r. This reduces the area to r^A. It also decreases the distance between the centers 
of any two disks by a factor of r. The second step is to translate each scaled disk back to its original position. 
During the translation the distance between the centers of any two disks increases continuously. By a result 
of BoUobas [6] , the union area of a set of congruent disks moving apart from one another continuously cannot 
decrease, and therefore the area of the final arrangement of scaled disks is at least r^A. □ 

5 Algorithm 

In this section we give an 0{n^) time algorithm to select and 3-colour a subset C of a set D of unit-disks 
so that the area bound given in Theorem [5] is realized. The proof of our theorem is constructive, and the 
only algorithmic issue is positioning the lattice so that at least lattice points are in Ul? (the union of 
all disks in D). We give an 0{n^) time algorithm for this lattice positioning problem. We also prove that 
the lattice positioning problem is 3SUM-hard, providing evidence that an 0{n?) time algorithm is the best 
we can expect with this approach. 

To solve the lattice positioning problem we follow the constructive approach used to prove Lemma [3l 
initially positioning the lattice arbitrarily and then translating all triangles in the lattice along with portions 
of UD to a fundamental cell F . Thus the one remaining detail is how to capture the translated portions 
of yjD so that we can compute a point of maximum depth. Our basic idea involves translating all of the 
disks and then computing and traversing their arrangement. Each disk d intersects at most 4 translates of 
F. We make 4 translated copies of d, and record which translate of F they come from. Computing this set 
of translated disks D' takes 0{n) time. Computing the arrangement of D', A{D'), takes 0{n?) time using 
the incremental insertion algorithm of Chazelle and Lee [8]. 

It is easy to traverse A{D') to compute maximum depth in D' — the depth increases when we enter a disk 
and decreases when we exit. However, this is not quite what we want; we want depth with respect to UD 
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translated to F, which is different from depth in D' due to disks that overlap originally in D. Our solution 
is to traverse A{D') maintaining the depth c,; in each translate i of F. Note that there are 0{n) translates 
of F that intersect UD. We also maintain a count c of the number of non-zero c^'s. The maximum value of 
c over cells of A{D') gives us what we want. 

We now prove that the lattice positioning problem discussed above is 3SUM-hard. A problem is 3SUM- 
hard if it is harder than the problem of determining whether a set 5 of n integers contains three elements 
a, 6, c € S such that a + b + c = 0. The best known algorithms for this problem take O(n^) and it is an open 
problem to do better [T3]. 

Theorem 6. The following problem is SSUM-hard: Given an integer k, real number s, and a set D of n 
unit-disks in the plane, determine whether a triangular lattice of side length s can be positioned such that it 
intersects the union area of D in at least k points. 

Proof. We show that our problem is harder than the known SSUM-hard problem of determining whether 
there is a point of depth fc in a set of unit radius disks in the plane. The more general problem for variable 
radius disks is proved 3SUM-hard in [4] and the reduction is easily modified to produce unit radius disks. 

Our reduction is as follows. Given a set D of unit radius disks in the plane, place an equilateral triangle 
T large enough to contain all of D. Expand T to a triangular lattice, and translate each disk of I? to a 
different cell in the lattice, with the same orientation. Let the translated set of disks be D' . Then there is a 
point of depth fc in if and only if the lattice can be translated to intersect D' in k points. This reduction 
takes linear time. □ 



6 Preliminary Bounds for /c-Colours, A: 7^ 3, and Future Work 

While we have focused on the case of 3 colours, since that is relevant for channel assignment in wireless 
networks, it is interesting to see what bounds can be derived for other numbers of colours. For 2 colours, we 
can use the techniques used to prove Theorems 1 and 2 but with a 2-coloured square lattice of side length 
2^/2. This approach leads to a bound of C2 ^ analogous to the bound from Theorem [1] and a weaker 
bound of C2 ~ 3^ for which an O(n^) algorithm exists for selecting disks to satisfy the bound. Full details 
are included in the thesis of Benjamin Lafreniere [TS]. 

For the more general problem of fc-colours. Theorem [7] presents some preliminary results, demonstrating 
a bound for all k such that k G {i'^ + ij + p \ i,j & N}. The number of such k lower than a given a; € N is 
given by Q[-^==), so the set of such k is thin (density 0). 

Theorem 7. Given k colours, where k G {i'^+ij+j'^ \ i,j ^ N} we can select and colour a subset of disks such 
that same-coloured disks are disjoint and their union area covers at least Ajj^^j^ where 6k = 

Proof. For all k € {i^ -I- ij -\- | i,j <E N}, distance Vk occurs in the unit triangular lattice L, and by 
^ rotational symmetry, an entire sublattice with side length ^/k exists. Thus we can partition L into k 
triangular lattices of side length ^/k and assign each a unique colour. We then scale L such that distance 2 
separates the enclosing disks of Voronoi cells of same-coloured lattice points by applying a scaling factor of 
afe — rr^ 9 ■ Now, each disk in D is assigned to the Voronoi cell containing its center. We select from each 

Voronoi cell one associated disk (if there are any) and colour it to match the Voronoi cell's lattice point. 
Note that by our scaling, same-coloured selected disks cannot intersect. 

If a point p is in [JD but is not in any selected disk, then the disk covering p intersects another disk with 
center in the same Voronoi cell, and the distance between their center points is less than the diameter of the 
Voronoi cell 5k ^ ^'^k- Now, if all selected disks were blown up by a factor of 1 -\- Sk, p would be covered 

by some selected disk and the union of selected disks would cover at most A{1 + Sk)'^- Thus, if we allow k 
colours, we can cover at least A ^^_^g^y2 . □ 
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A Proof of Lemma [H 



Our first claim is that the minimum area of intersection is achieved by a disk d positioned such that X, 
the center point of hexagon h, lies on its boundary. Suppose this is not the case. By symmetry, it suffices 
to consider possible placements of Y, the center point of d, within the intersection of wedge BXK and the 
unit-disk centered at X in Figure [SI For any position of Y , moving Y to the right along a line parallel to AB 
decreases the area of intersection, since the portion oi d — h lying above the supporting line of AB stays the 
same, and the portion oi d ~ h below the supporting line of AB strictly increases (by containment). Thus 
we can move Y to the right until it lies either on XK or the boundary of the disk centered at X. For Y on 
XK , moving Y toward K decreases d n /i because when the diameter of d parallel to BC lies strictly inside 
h, the area oi d — h increases (by containment), and when the diameter is not strictly contained in h, the 
area of dn /i decreases (by containment). 




Figure 9: Sliding disk d outward from the center of h results in a smaller intersection between d and h 

Thus we can restrict our attention to the minimum area of intersection with h among disks whose 
boundary contains point X. To find this minimum, we assume that X — (0, 0) and express the area of 
intersection f{0) in terms of angle 9 between the center of a disk d, the center X of /i, and the x-axis. There 
are two general cases to consider. Case 1 occurs when d contains two vertices of h (e.g. Figure [TU)) . Case 2 
occurs when d only contains a single vertex of h (e.g. Figure [TT|) . Note that in either case we can express 
the area of intersection as the sum of the area of a polygon and a circle sector. For instance, in Figure [TO] 
the area of intersection is the sum of the area of polygon ABCED and the area of the sector of d interior to 
angle ABC. 

By symmetry, we need only consider the area of intersection for < 9 < ^ . As a result, we can 
use the cases shown in Figures ITOl and ITT] to derive a formula for f{9). Specifically, we use the symbolic 
geometry package Geometry Expressions to derive formulas relating 9 and the intersection points between 
the boundaries of d and h (i.e. points A and C in Figures [TOl and \TT^ . The derived formulas for these points, 
along with formulas for the other points in Figures [TOl and [TT] are given in Appendix I A . 1 1 and I A . 2 1 respectivelv. 
From these formulas we express the area of intersection in terms of 9 using standard formulas for the area of 
polygons and circle sectors. This gives us formula fi{9) for the area of intersection for < 9 < arcc os(|) — ^ 
and formula f2{9) for the area of intersection for arccos(|) — ^ < 9 < ^ (see Appendix I A. II and I A. 21 By 
symmetry, we extend this to formula f{9) for < < ^ given in Appendix lA. 31 

A plot showing f{d) for the interval < 9 < ^ is given in Figure [T2l Using Maple we find that f'{6) (the 
first derivative of f{9)) is for 9 = ^. Thus, by symmetry, the minimum intersection occurs when X, the 
center point of d and a vertex of of h are collinear. Computing the value of f{9) at any one of these points 
gives our value for A, specifically 
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Figure 10: Calculating the intersection area as a function of angle 9 (Case 1). 
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Figure 12: Plot showing the area of intersection for < 6* < ^. 
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A.l Derived Formula for < 6 < arccos(|) — | (see Fig. [TOl) 

1 l//2 1 1 3 

- - - J 1 - ( - - + - sin (6») + - cos (e) j - - sin (0) + - cos (0) 

l//2 1 1 X^^l 1 
- W 1 - I - - V3 + - sin (6») + - cos (e) ) + - sin (0) cos (6*) V3 



■■ cos(6») 
: sin{e) 



1 l//2 1 1 3 
-v^ Wl-(-v^+ - sin {9)V3 cos (6») J + - sin (9) Vs + - cos (9) 

32yy32 2 y 4 4 

1//2 1 1 1 
- W 1 - ( - \/3 + - sin (6*) \/3 cos (6*) ) + - sin (6) + - cos (9) 

2yy3 2 2 y4 4 



3 

-1 - 



V3 



V3 

2 

3 



= - (A4e){Byie) - Dy(e)) + B^(e){Cy{e) ~ Ay{e)) + C4e)(i?^(e) - By{e)) + E^{e){Dy{e) - Cy{e)) + D4e){Ay{e) - Ey{e))) 



1 / ^^^^^^ [ -{~B.{e) + c^e)) (-A^(e) + Byje)) + ~ Cy(e)) (A^je) - B^e)) 

2 ''''' ""^ V + c,(0)) - B,(0)) + {By{e) - Cy(e)) {-Ay{e) + By(e)) 

A. 2 Derived Formula for arccos(|) — f < < f (see Fig. [TTl) 



)) 



By(e 
1 



- W 1 - I - - + - sin (6») + - cos (e) ) sin (e) + - cos (6») 

32y\,32 2 J 4 4 



1 



2 

: C0S(6») 

: sin{0) 

2 r- 
■ - V3 



^ ^1 - (- ^ + ^ sin (e) + i cos (e)) + i sin (0) - 1 cos (5) 



- i ^-3 + 12 cos (e) - 9 (cos (6»))2 + sin (5) 



;V3 



f2(e) = - (A4e)(Bj,(0) - Dy{e)) + B^(e)(Cy{e) - Ay(e)) + c4e)iDy{e) - By(e)) + D^{e)(Ay(e) - Cy(e))) 

1 / ^^^^^^ ( ~{-B.{e) + c^(e)) {-Ay(e) + By(e)) + (B^e) - c^e)) (A^e) - B^e)) 
2^ V {-B^(e) + c^:{e))(A^{e)-B^{e)) + {By{e)-Cy{e)){-Ay(e) + By{e)) 



A. 3 Derived Formula for O < 9 < ^ 



fid) 




if < 9 < arccos(|) - | 
if arccos(|) - f < 9 < | 



if arccos(^) 
if "I < 6» < arccos(| 
if arccos(^) < 6 < 



15 



